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. The recent optical transmission experiment of Karrai et al. has probed the local- 
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Q\ , paper we develop a microscopic description of single vortex dynamics, based on the 



ized quasiparticle excitations of vortices in the superconductor YBa2Cu307. In this 



Bogoliubov-deGennes equations and self-consistency through the gap equation, to 
determine the response of vortices to a transverse time dependent electric field. It is 
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applicable to the low temperature, clean, extreme type II limit. In the limit of large 

O 

planar mass anisotropy it simplifies. Thus it may be especially relevent to materials 
such as NbSe2 and high temperature superconductors. Of particular interest is the 
response of the vortex at frequencies near the minigap, A 2 /Ep, where A is the bulk 
energy gap and Ep is the fermi energy. A dissipative equation of motion for vortex 
cores valid at non-zero frequencies is derived. We give a clear microscopic meaning 
to the vortex drag parameter. The expected dipole transition between quasi-particle 
states localized at the core is hidden because of the self-consistent nature of the vor- 
tex potential. Instead the vortex itself moves and has a resonance at the frequency of 
the transition. We calculate the conductivity of vortices as a function of frequency. 
A analogy is made to the Mattis-Bardeen result for the electrodynamic response of 
bulk superconductors that, unless translation invariance is broken, single particle 



properties are 'invisible' to a long-wavelength probe. However we show that upon 
adding a translation invariance breaking term the dipole transition re-appears. This 
approach may eventually form the basis of a microscopic theory of vortex pinning at 
non-zero frequency. 
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I. INTRODUCTION 



The condensed state of type II superconductors can possess a topological defect, a vortex, 
at which there is a zero in the order parameter. The existence of these defects is made 
energetically favourable upon the application of a sufficiently large magnetic field or current. 
Associated with these defects are quasiparticle states with wavefunctions localized near the 
vortex cores and energies within the bulk energy gap. These states are superpositions of 
electrons and holes and one may think of them as trapped in the vortex by continual Andreev 
reflections due to the spatially varying order parameter. Their quantitative description in 
a simple limit was formulated some time ago. Caroli, deGennes and Matricon |lj and later 
Bardeen et al. [g] calculated the energies and wavefunctions of these discrete levels using the 
Bogoliubov-deGennes (BdG) equations. 

The electrodynamical response of type II superconductors in the mixed state is deter- 
mined by the dynamics of vortices. It is usual to approach this by saying that the phase slip 
from vortex motion is associated with voltage drops and hence gives rise to dissipation in the 
superconducting state. Another, less conventional, way to look at things is to note that the 
low energy excitations of the system are associated with the states inside the bulk energy 
gap which are localized near vortex cores. For most practical applications, however, vortex 
pinning is the most important effect. But in the case of very clean systems (possibly probed 
at high frequencies) one is interested in dissipative flux flow. The microscopic calculations 
mentioned above helped form a basis for theories of vortex motion based on the idea of a 
'normal' core such as those of Bardeen and Stephen |§, Nozieres and Vinen [Q and others 

i- 

Even in the presence of weak pinning, flux flow may be relevant if pinning is overcome 
at high enough temperatures (but below T c ) or high enough frequencies. Gittleman and 
Rosenblum ||[7j studied dissipation in the mixed state of conventional superconductors as 
a function of frequency and noticed that there was a crossover to essentially free vortex 
behaviour over one or two orders of magnitude centered about some frequency. What is 
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significant for this work is that the crossover frequency region appears to be of the same 
order of magnitude as the expected spacing between discrete levels in vortex cores. That 
spacing is of order A 2 /Ep or ft 2 /m£ 2 (where A is the bulk gap, Ep is the Fermi energy and £ is 
the coherence length) and follows from applying the uncertainty principle to a quasiparticle 
confined to a vortex core of size £. Unless impurities radically alter the vortex size this 
should be the characteristic energy scale. The comparison is made in table H where some 
unpublished data of Schleger and Hardy || are included. The fact that the de-pinning and 
level spacing frequencies are roughly comparable means that any microscopic understanding 
of the de-pinning transition as a function of frequency should take into account the discrete 
spectrum of states inside vortex cores. 

Caroli and Matricon || discussed possible ways of observing the discrete structure within 
vortex cores through ultrasonic and nuclear magnetic relaxation. None of these methods has 
so far succeeded. A few years ago Hess et al. probed the gross structure of vortices in 
NbSe2 by scanneling tunneling microscopy (STM). The results were qualitatively explained 
by theories based upon the Eilenberger equation approach of Kramer and Pesch |12|| . The 



STM experiments however had a resolution of about O.lmeV whereas the energy separation 
of the discrete levels is expected to be about lOmK. STM experiments probing the vortex 
cores of high temperature superconductors (where the energy scale is expected to be of 
order lmeV or more) would be extremely desirable but have not yet been accomplished. In 
passing we note that in charge density wave systems similar microscopic vortex structure 
and possible experimental implications were proposed by Maki and Huang ||10|| . 

Recent experiments carried out by Karrai et al. |14| have measured the optical trans- 



mission of the mixed state of YBa2Cus07 at low temperature. They observed a knee in the 
normalized transmission coefficient in the frequency range 50 — 100cm" 1 and attributed it to 
a quasiparticle resonance. Some features are of interest. First the resonance frequency is at 
least 3 or 4 times higher than expected from a simple microscopic calculation. Second, the 
fitted relaxation rate of the resonance is very large: of the same order of magnitude as the 
resonant frequency itself, 50cm 1 ~ 10 12 s 1 . The is roughly the same order of magni- 



tude as the relaxation rate of YBa2Cu30y as measured by microwave surface impedance [15 



The experiment is performed at a temperature T = 2K, over an order of magnitude lower 
than the resonant frequency which allows us to deal with only the lowest energy excitations. 
One would like to know whether the observed effect is due to the intrinsic response of free 
vortices, or whether it is more important to consider defects, with which the vortex cores 
would be highly correlated. 

In a previous paper |16| the response of quasiparticles in cores of vortices (in the clean, 
low temperature limit) to a long wavelength electromagnetic field was considered. It was 
concluded that the motion of the vortex itself had to be taken into account because it 
was a self-consistent potential. In this paper this idea is re-derived in a different fashion 
and expanded in order to be directly relevant to infrared measurements in high temperature 
superconductors. The paper is organized as follows: In section [H| we introduce the formalism 



associated with the Bogoliubov quasiparticles and the BdG equation. In sections [III] and 
|TV| we describe displaced and moving vortices respectively using this formalism. In section 
[V] we derive the equation of motion, calculate the response of vortices to a long wavelength 
electromagnetic wave and discuss its consequences. We discuss some effects of pinning in 
section [VTL It is followed by a discussion and summary in sections |V11| and |V111|. 



II. MICROSCOPIC FORMALISM 

The full Eliashberg-Gorkov formalism of BCS theory is, of course, intractable for the 
problem of a moving vortex. Instead we consider the simpler BdG equation for s-wave 
superconductors. This equation is valid in the case of a short range instantaneous pairing 
interaction. It may also be derived in the quasiclassical limit kp^ » 1 where is the 
fermi momentum and £ is the coherence length. That limit is the leading order result in 
the formulation of Eilenberger JET]]. We shall supplement the BdG equations with a local 
gap equation A (r) = V(c-f (r) C| (r)) which determines how the vortex itself moves (self- 
consistently). Of course high temperature superconductors stretch the assumptions made 
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here. The coherence length is only about 5-10 lattice constants parallel to the Cu-0 planes 
so that kpC, >> 1 is barely satisfied. The discreteness of the lattice may also come into 
play. Real superconductors are of course non-local and high temperature superconductors 
may even have d-type symmetry of the order parameter (for which the microscopic structure 
of a vortex has not yet been calculated). Also, in general, the pairing interaction may be 
retarded but since we are concerned with very low frequencies relative to the gap, some 
relevence to real materials may remain. 

The eigenf unctions for Bogoliubov quasiparticles are two component objects which can 
be thought of as the electron and hole amplitudes of the quasiparticle, 
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(2.1) 



ip satisfies a Schrodinger equation which is just the BdG equation, 
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where a z is a Pauli matrix, Ep is the chemical potential, and A = |A(r — 
r )| exp (— i9(r — r )). 9(r — r ) is the angle about the center of the vortex r measured 
from the x axis. By convention the z direction shall be in the direction parallel to the 
vortex. 

Going to second quantized formalism the conventional definition || of the quasiparticle 
creation operators is, 
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(2.3) 



cj(r) -c T (rj j 

where \x is an index for the low energy solutions of this equation and Cf, denote spin 
up and spin down electron operators respectively. It is conventional to define the creation 
operators only for positive energy states (they are a complete set of states). However it 
will be convenient for us to use the complete set of states which includes both positive and 
negative energy states but eliminates the spin degeneracy. That is, we shall use, for > 0, 
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7^rt ? , (2.4) 

and for e M < 0, 

7M = ri M , (2.5) 

where —\i refers to the corresponding time reversed positive energy state. This is possible 
because if (u^ v^) is a solution for energy then (v* — «*) is a solution for energy — e^. 
Finally, as we shall see below, we have a ladder of states from negative to positive energy 
which are filled for e M < 0. The inverse transformation from electrons to quasiparticles is 
then 

4(r)=E M 7K( r ) > 

cftr) =E,7W r ) • ( 2 - 6 ) 

We shall be interested in the extreme type II limit with H « H c2 . This can certainly 
be satisfied in the high temperature superconductors. Vortices are well separated and there 
is no quasiparticle tunneling between vortices (although see reference |2I| for a treatment 
of that situation). The magnetic field applied to create the vortices may be ignored in the 
BdG equation. Because the magnetic field is spread out over an area A 2 it's importance 
compared to the phase of the order parameter is reduced by £ 2 /A 2 where £ is the coherence 
length and A is the penetration depth. 

The low energy eigenf unctions for fixed k z , \i « kp±£, and the radial coordinate r << £ 



arc 



^2 



(2.7) 



with angular momentum index \i = ±|, ±|, ... and kp± referring to the Fermi momentum 
projected onto the k x , k y plane. At distances of order £ from the vortex center the wavefunc- 
tions begin to decay exponentially Jl]U. The energies as calculated by Kramer and Pesch 
T2] |, who analytically accounted for some self-consistency effects due to the gap equation 



A(r) = \/(c T (r)q(r)), are 



s = ^ ln(* 
/cfVfcos 2 9 2 



ln(-£„cose/fi), cosG = k F± /k F . (2.8) 



The zero of energy is the Fermi energy. 

The logarithmic factor is a simple rescaling of energies. It was derived by Kramer and 
Pesch and arises from self-consistency of the vortex solution. At low temperature, energy 
is favoured over entropy and so the vortex pinches together. This increases level spacing, 
pushing down the energies of occupied states and reducing the occupation of excited states. 
This factor will not be important for this paper but its effects were explored by Bardeen 



and Sherman ]18| and Larkin and Ovchinnikov 

In this paper, in order to simplify matters, we will suppose that the Fermi surface is open 
and nearly cylindrical. That is, the system is nearly two dimensional and the dispersion in 
the (by convention z) direction perpendicular to the planes is very small, i.e. the vortex 
is aligned perpendicular to the planes. This is a rather good approximation for materials 
such as NbSe2 and the copper-oxide high temperature superconductors. Moreover there 
is another factor which helps. Near two- dimensionality is equivalent to k F ± = kpcosQ 
varying little with k z . As one can see from Eq. (|2.8|) , the dispersion in the k z direction 
enters through the cos0 in the denominator. Each angular momentum level broadens into 
a band. Nevertheless, the one- dimensional density of states associated with the variation of 
this cosine diverges at O ~ and this strengthens the approximation of neglecting the k z 
dispersion. 



III. DISPLACED VORTEX 

The quasiparticle states are a complete set of states so that once the density matrix in 
that basis is specified the state of the system is defined. Conversely it is possible to describe 
a displaced vortex in this basis because the order parameter is completely specified in terms 
of the quasiparticles through the gap equation. By using Eqn. ( |2.6|) and substituting into 
the gap equation we have, 



S 



As an aside recall that in the well-known Ginzburg-Landau solution for the order-parameter 
profile of a vortex core in an s-wave superconductor we have A(r) ~ r, as r — > 0. From the 
series representation of the Bessel functions, J n (x) ~ x'™' for small x and integer n. From 
Eqs. ( p.7|) and (|3.1j) we can see that this property is satisfied by the vortex on a microscopic 
level, but it is only the /i = ±1/2 states that contribute to the linear in r component of the 
order parameter. At temperatures higher than A 2 jksEp the occupation of fi = ±1/2 states 
is about the same and the linear in r component is washed out because J- n (x) = (—l) n J n (x). 

Suppose the vortex is displaced by a small amount <5r in a direction which makes angle 
0o with the x-axis. Let the occupation be diagonal before displacement, that is, (7^7^) = 
S tlv f(e L ,), where /(e) is the fermi distribution function. The changes Su^r) = — <5r ■ Vw^r) 
(and similarly for f^(r)) may be approximated by assuming that u and v behave roughly like 
Bessel functions. This is explained in detail in the appendix. In the undisplaced coordinates, 
to linear order in 5ro, and for small r, 

SA(t) « - V s -^ £„ (f(e u ) - f(e u+1 )) 



x 



e ifr v * ( r )«,(r)±e-^°<(r)w, +1 (r) . (3.2) 



It is important to note that small displacements can be described in terms of changes in 
the quasiparticle occupation or density matrix. The above displacement can be represented 
by the change 

^ +1 7,) = ^^[/(^)-/(6, + i)]e^, Vu (3.3) 

(and the Hermitian conjugate) in the quasiparticle density matrix using the undisplaced 
basis. 



IV. MOVING VORTEX 



In this section we consider a vortex moving with a small velocity and show, using a 
Galilean transformation, how it can be approximately described by quasiparticle excitations 
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in the stationary frame. Suppose that the vortex moves with velocity v. Then using the 
invariance of the BdG equation the eigenf unctions transform like 

u(r) -> u(r - v£)e imv ' r 
v(r) v(r - vt)e- imv - r . (4.1) 



The gap transforms like 



A(r) -> A(r - vt)e 2mv r , (4.2) 



while the chemical potential transforms as 

1 e 

E F -> E F + -mv 2 - -v • A , (4.3) 

2 c 

assuming the gauge field A is constant. From the usual gauge invariance there is an arbi- 
trariness as to whether to put the extra phase in A or A. For what follows we shall take 
A = 0. 

We wish to evaluate a general component {'y v j r i r ) v ) of the density matrix in the stationary 
frame. But the quasiparticle operators (for which the corresponding density matrix is taken 
to be diagonal) in the moving frame look like 

/yl = j d 3 rc}(rK(r - vt)e imv r + q(r)^(r - vt) e - fmv - r (4.4) 

for v > and 

fyl = J d 3 r - c\{r)v*_ u (r - vt)e imv - r + Cl (r)u*_ u (r - wt)e- imv - r (4.5) 

for v < 0. Let us consider a vortex which is moving but not displaced. That is, we take 
t = in the above expressions. Also we retain only terms up to first order in the velocity. 
Using the definitions of Eqs. (|4.4 ) and (|4.5|) and expanding the exponential we can rewrite 



the anti-commutator as 

{7^} ~ V + / d 3 vipl(r)a z (imv ■ r)^(r) . (4.6) 
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Using the results of appendix, we may evaluate the matrix element in Eq. ( |4.6| ) and find 

where 6q is the angle that v makes with x. 

Now we are able to calculate the density matrix for a moving vortex, but in the basis of 
eigenstates of a stationary vortex. To linear order in the velocity we have 

(7h,±i> « T[f{e n ) ~ fK ±l )} V F±e . (4.8) 

4t(e u ±i — e u ) 



Note the similarity in form of this expression with that of a displaced vortex, Eq. fl3.3| ). 
Comparing the two we find that the corresponding displacement is 5r = hv/(A 2 / Ep) at 
an angle 7r/2 relative to the velocity direction. The naive gap equation is not satisfied for 
a moving vortex! The quasiparticles are displaced relative to the center of the vortex and 
have a non-equilibrium distribution as seen from the lattice frame. 

At this point we should note the effect of dispersion in the k z direction on the energy 
denominator in Eq. (O). From Eq. (2.8) we see that it varies as cos _2 as we change k z . 



Therefore the displacement discussed above is not the same for states of different k z . In fact 
it will be clear that the different k z components of the vortex will not move together. We 
shall ignore this for now and assume a rigid motion of the vortex but this is one place where 
a relaxation of the two-dimensionality assumption will complicate matters. Vortex bending 
is another matter and will be discussed in section |V| in conjunction with pinning. 

V. RESPONSE TO EXTERNAL FIELD 
A. VORTEX EQUATION OF MOTION 

We shall be interested in the linear response of a vortex to a uniform time- varying super- 



current or electromagnetic wave. In a previous paper |16| an equation of motion was derived 
for a vortex in the presence of this perturbation. The derivation was a bit cumbersome 
because it relied on tracking the evolution of the quasiparticle density-matrix as a function 
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of time for small times and then extrapolating to all times. In this paper we shall derive 
the same result by an entirely different means which takes advantage of the results derived 
in the previous two sections. 

Consider a long wavelength electromagnetic wave incident in the z direction, A, with 
polarization making angle 6q with the x axis and perpendicular to the vortex line itself. The 
perturbation to the Hamiltonian is 

d 1 



e the [ 
— A • p = A < sin 

mc mc 



*-(*o + f) 



d 

TT" + COS 

or 



0-(0 o +f: 



rdd 



(5.1) 



The vortex will move with some velocity in the presence of a background superfluid 
velocity vg = -eA/mc resulting from the applied field. These velocities will be assumed 
uniform in the z direction (along the length of the vortex). Such an assumption would be 
valid if the distance that the wave penetrates the superconductor (the smaller of the London 
penetration depth and the skin depth) were long and the distance between pins were also 
long compared to the coherence length. 

First let us present a simple example to illustrate how we shall calculate the motion of 
the vortex in this paper. The basic idea is to calculate time derivatives of the quasiparticle 
density-matrix and use the results of sections |TTT| and [TV] to identify the corresponding 
motions of the vortex. 

From the appendix we find that the matrix element of Eq. (|5.1|) is approximately 



/ rf 3 r< ±1 (-— A ■ pty, = £^i e ^o + f) = W (5 . 2) 
J p mc 2mc 

If A were time dependent and this were the only perturbation then we would have 

d -H%-x)(t) = [f(e») - f(e^)] x i-i^-) (5.3) 



and its Hermitian conjugate. 

Compare this result with Eq. ( |3.3j ) and assume that there is a time dependence in 8tq. 
Making the identification d(8ro)/dt = we see that = vg. This is the well known 
result, that, in the absence of dissipation or other forces the vortex moves with the same 
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velocity as the background superfluid. This is the expected result from Galilean invariance 
or from the fact that in wave mechanics, a uniform gauge field, A, boosts the group velocity 
of all waves by a velocity V5 = — (e/mc)A. 

Let us now proceed with the derivation of the full equation of motion. Suppose that the 
incident electromagnetic field A is time dependent but still uniform. It will be convenient 
to perform a gauge/ Galilean transformation which is a boost by velocity V5 and eliminates 
A. Under this transformation, 

u(r) -> u(j)e^ e / h ^ A - r 
v(t) -> v ( r )e- i( - e/hc)A - r , 

A(r) -> A(r)e 2i ^ hc ^ A - r . (5.4) 

The time dependence of A rewritten as S = A/c results in a perturbation a z eS ■ r. Now 
we are left with A = and a vortex which has velocity — V5. The quasiparticle density 
matrix has off-diagonal elements 

(Th^i) - =FL/(fr) - f(e v±1 )} |VL 9 ~ Vs|fejr tT° • ( 5 - 5 ) 

M\£v±i — ^v) 

The angle 0o is the angle of — V5 with respect to the x axis. Here we remark that 
the Fermion occupation f{e v ) does not change to linear order in velocity. Because of the 
presence of a gap the energy goes quadratically with the velocity. 

Now the time evolution of the off-diagonal elements of the density matrix coming from 
the difference e v 7^ e v ±\ corresponds to the motion of the vortex at velocity vx — V5 as 
derived in section [IV|. However a moving vortex also has a relative displacement between 
the order parameter and the quasiparticles relative to what the gap equation for a stationary 
vortex would give. This displacement is 

Sr = (5.6) 

in a perpendicular direction 0' o = <po — vr/2. Thus the displaced vortex itself will affect the 
quasiparticle density matrix through a perturbation 5 A = —5r' Q ■ VA. Its matrix element is 
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w, 



III' 



d 3 r 



1 5A^ 



V 



SA* 

In the appendix we show that it has the value 

Sr' Q k F ± 



(5.7) 



W, 



2i 



(5.8) 



The contribution of the relative displacement <5i"q to the evolution of the density matrix 



is 



|(t!7m-i> = Lfe) " x {-^-^} = ^'|vl - v s | k F± e**>. (5.9) 

Comparing with Eq. ( |4.8| ) which relates velocity to off-diagonal elements of the density 
matrix we may interpret the results in terms of an acceleration 



a = |v L - v s |(e„ + i - e„)//l 



(5.10) 



at an angle <po — tc/2, or, in other words, 

A 2 



hE f 



(Vi-V s )XZ . 



(5.11) 



As for the electric field we require the matrix element (using the results from the ap- 
pendix) 



ijAe£ ■ ra z ijj v 



eSk 



2m(e u+1 - e v ) 

where <po is now the angle of S. From the electric field then we have a contribution 



(5.12) 



d 



eSk 



1 e^ . 



(5.13) 



The corresponding acceleration, by comparing with Eq. (|4.8|) is 



eS . 
a = = v 5 

m 



(5.14) 



The last influence on the vortex motion we shall consider in detail is dissipation. Suppose 
that we are in the lattice rest frame and looking at vortex moving with velocity v^. We 
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see that there is apparently a non-equilibrium off-diagonal component whose magnitude is 
proportional to Vl and whose phase gives the direction of w L . Let us assume that, due to 
scattering from things that break translation invariance such as impurities or phonons, this 
off-diagonal component has a lifetime r over which it decays and re-appears as a thermalized 
(diagonal), isotropic contribution to the density matrix (i.e. not carrying current). In other 
words we assume that the core can maintain thermal equilibrium with the lattice unlike, 
say, the case of 3 He in which the cores of vortices can heat up. The parameter r which will 
appear in the vortex equation of motion in the drag term and in the conductivity has a clear 
microscopic meaning. It is the lifetime of low energy quasiparticle states. 

This theory of the dissipation then gives an additional component to the time derivative 
of the vortex velocity, 

v L = --v L . (5.15) 

T 

Putting all the contributions together we have the final equation of motion, 

A 2 1 
vl = v s + — ^-(v L - v s )xl v L . (5.16) 

It is useful to compare the present expression with that of Nozieres-Vinen ||. Their 
equation is derived by balancing forces. There is the 'magnus' force (/m/2)(vs — vl) x z 
where n is the (superfluid) number density of electrons. Comparing this to the corresponding 
coefficient in the present equation of motion, taking n = k^ ± /2n and the in-plane coherence 
length = hvp^/nA, one may extract a 'mass' of the vortex M ~ m±(kF±l;±) 2 /4 per 
unit length in the z direction. This expression is perhaps a microscopic justification for the 
concept of a 'normal core' of size £, even at very low temperatures when there is a gap in 
the density of states of single particles. This is should be contrasted with another definition 



of the mass of a vortex as recently discussed by Duan and Legget pifl . This definition is 
based upon the change in energy of a moving vortex at order v 2 which has been ignored in 
this paper. 

Let us pause and summarize the situation. We began by asking whether the applica- 
tion of a long-wavelength electromagnetic field could cause dipole transitions between the 
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discrete quasiparticle states in the cores of vortices. Bearing in mind the approximation of 
rigid vortex motion and other assumptions we find that an applied electric field, instead of 
causing dipole transitions, causes the density matrix to evolve off-diagonal matrix elements 
corresponding to vortex motion itself (after applying the gap equation). The vortex does 
not stand still and allow a dipole transition to take place, as the core of an atom. The 
difference is that the vortex is a self- consistent potential. Thus the experimental observation 
of such quasiparticle resonances will not be a simple dipole resonance (except in the pres- 
ence of pins which will be discussed in section [VI]). The STM tunneling experiment differs 
from the present case in that the tunneling process introduces an extra particle which is not 
correlated with the quasiparticles in the vortex whilst the effect of an applied electric field 
is to change the off-diagonal components of the density matrix. 

A simple understanding of the equation of motion comes from looking at the homogeneous 
solution. It can be written 



VLx 

\ VLy ) 



-t/T 



a + e 



iQot 



(A 



-iQot 



( I \ 



v -v 



(5.17) 



which corresponds to moving in circles, with a definite sense of rotation, at a frequency 
Qq = A 2 /Ep- This is the sort of state one hopes to excite with the external probe. Therefore 
experimentally it would be important to probe this system with polarized waves In the next 
subsection we discuss what happens when we do so. 



B. DISSIPATION IN A SINGLE VORTEX 

It is very illuminating to calculate the dissipation in a single vortex as a function of 
frequency and polarization for constant magnitude vs- In the next subsection we shall 
calculate the conductivity for a finite density of vortices. First it is necessary to extract the 
paramagnetic current carried by the vortex core. The paramagnetic current is 

J(r) = E^W t c,(r) . (5.18) 
16 



Substituting for the electron operators using the inverse transformation Eq. (|2.6|) , and using 
the results of the appendix to evaluate matrix elements of V we find that 

/ ft • J = ^ E [lUil^* - 7 W°) , (5-19) 

where n is some arbitrary unit vector pointing in direction <p . 

Now we may calculate the (spatially averaged or zero momentum) paramagnetic current 
associated with a vortex motion by comparing Eqs. (|5.19|) and fl4.8Q . It is 

per unit length in the z direction (for a non-cylindrical Fermi surface we would average over 
k z ). 

Let us calculate now the dissipation as a function of frequency for a single vortex. For 
simplicity we shall keep the amplitude v s constant but we shall allow it to have an arbitrary 
polarization. The steady state solution of Eq. ( p,16| ) for v s (t) = V5(0)exp {iut) is 

Qqtvsx - (1 + iur)v Sy 



VLy - V S y 



(5.21) 



(1 + turf + (ft r) 2 
and another equation with x,y interchanged and —* —flo- 

Eq. ( |5.16| ) at zero frequency was introduced by deGennes and Matricon |22|. This had 



the drawback of not allowing for small conductivities seen in flux flow experiments [^] and 
was discarded by Nozieres and Vinen in favour of an equation of motion where the dissipation 
acts on vj rather than v^. Nevertheless the solution Eq. ( |5.21|) agrees with the Nozieres- 
Vinen equation j|] at low frequency and low dissipation. It is probably not valid except in 
the clean limit where the levels are well defined. Moreover, because of the the frequency 
dependent term there is a 'resonance' at uj ~ Qq with a width r _1 . If one simply added 
a vs term to the Nozieres- Vinen equation one would obtain an unphysical divergence at 
that frequency. In addition, even though Eq. ( p,16| ) was not derived in the large dissipation 



limit, that limit makes sense for the present equation of motion. The vortex simply stops 
moving. The Nozieres- Vinen equation and its relatives have divergent behaviour because in 
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these models the current through the vortex core is forced to be equal to the background 
transport velocity V5. 

There are two contributions to the dissipation which come from a current (either 
the background superfluid or that due to vortex motion) that is in phase with an elec- 
tric field (time derivative of the supercurrent or transverse motion of the vortex). The 
first is transverse vortex motion in phase with the supercurrent. For clarity let v$ y = 
vsx exp 16 = vs/ V2- The parameter 9 simply controls the polarization. There is an in- 
duced voltage per vortex z x v L (h/2e). The supercurrent density vsne gives dissipation 
— N v (hn/2)Re(vl x vs y — v* Ly Vs x )i where N v is the vortex density. Using Eq. (|5.21|) 



Mv Lx VSy ~ V Ly V Sx ) = -«§»0T [1 _ K)a V ( / >QT } a] aH ^ ((riT)a 

l + (^r) 2 -ffi r) 2 
[l-(o;r) 2 +(nor) 2 ] 2 +4(wr) 2 ' 



I "jn/??; 2 ,-.,-r l+(^) 2 -ffi r) 2 (f, 99) 



This first term contributes at low frequencies (relative to A 2 /Ep) and corresponds to the 
usual phase slip dissipation mechanism for flux flow. 

The second source is the current due to vortex motion which is in phase and parallel with 
the applied electric field. A straightforward calculation gives the average current density due 
to vortex motion to be 2vi{Ep / 'A) 2 N v e per unit length in the z direction. With the electric 
field £ = {m/e)vs the dissipation is 2m(E F / A) 2 N v Re(vl y iuvsy + v* Lx ioJv Sx )- Using Eq. 
flOU ), we obtain |4 



r> i ■ i 2 2 1+(u)t)' 2 — (Slnr) 2 

Re{v Ly zujvsy + v Lx iuv Sx ) = v 2 s uj 2 t [1 _^ + ( noT )2^ {u]T) 2 



-sm9v 2 s Q ujT [1 _ ( J T) 2+(n t ( )"]2+ ) 4 (ci;T) 2 • (5-23) 

This second part contributes mostly to high frequencies. It may be understood better 
by letting Q — > 0. In that case we recover a Drude-like expression 1/[1 + (^t) 2 ] multiplied 
by the volume of the vortex cores. This analogous to the well known zero frequency result 
for flux flow which is that as H — > H C 2, the resistance in the vortex cores multiplied by the 
core volume matches the resistance for the normal state above H C 2. 
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For clarity and in order to add together these two terms we shall make the oversimplified 
but well-defined choice that superfluid density be equal to the density of all electrons, n = 
k F± /2% (again, per unit length in the z direction) and Ep = h 2 k F± /2m independent of k z . 

With this definite choice we may then plot the dissipation. In Fig. (|1|) we have plot- 
ted the dissipation for right and left circularly polarized V5. For one of the polarizations 
we see dissipation but no special behaviour near the characteristic frequency of about Qq 
(depending on the magnitude of r). In the other polarization the vortex does respond at 
the characteristic frequency but essentially as an anti-resonance. There is a minimum in 
dissipation. 

C. CONDUCTIVITY 

In this subsection we calculate the experimentally accessible conductivity at momentum 
q = 0. The usual straightforward way of calculating the conductivity would be to take the 
gauge invariant current-current correlation operator, re-express it in terms of Bogoliubov 
operators using Eq. ( |2.6| ) (allowing for time dependence of the fermion operators) and 
evaluate its expectation value. As usual the (time) Fourier transform of this object will 
have poles at the excitation frequencies of the system. The lowest quasi-particle-quasi-hole 
excitation would contribute a pole ostensibly around Qq. That would lead to a contribution 
proportional to 5{uj — Qq) in the real part of the conductivity. 

However that is not the whole story because of the 'residual' or 'final-state' interaction 
between these quasiparticles. It can modify the energy from what one expects given the 
effective single-quasiparticle energy spectrum. Very generally, when symmetry breaking and 
long-range order occurs in a Fermionic system resulting in an energy gap for single-particle 
excitations these residual interactions can mix particle-hole excitations with the collective 
(Goldstone) mode by creating a bound state in the particle-hole channel. Another way to 
look at this is to say that while a mean field or Hartree order parameter is constructed so 
that single particle excitations have positive energy and do not scatter, when a particle and 
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a hole are present the change in the order parameter due to the hole affects the particle and 
vice versa. 

Since the position of a vortex is arbitrary (the real magnetic field that creates it is arbi- 
trarily uniform), the overall energy of a vortex should be unchanged if it is translated by some 
amount. If this amount is small compared to the coherence length (and this is certainly the 
case in the experimental situation), then it can be described by one quasiparticle-quasihole 
pair as given in section |HI[ Thus one expects from translational invariance that the residual 
interaction will reduce the quasiparticle-quasihole excitation energy to nearly zero. 

Unfortunately, an explicit demonstration of final-state effects is difficult because in order 
to get precisely zero energy one must possess precise self-consistent order parameter and 
single-particle states. This calculation would require detailed numerical work even for a 
simple zero-range pairing interaction and would be hopeless for any realistic model. Nev- 
ertheless, let us show, analytically, in a simple model, that the correction to the excitation 
energy is of order — A 2 /(i?pln[2u; c /A]). We begin with a simple zero-range pairing interac- 
tion 

- V'y'drc{(r)c ir (r)cJ(r)c 1 (r) . (5.24) 

In this case the residual interaction is 

-Vjdr (c{(r)4(r) - A*(r)) (q(r)c T (r) - A(r)) . (5.25) 

This expression, converted to Bogoliubov operators, is diagonal for the simplest // = — 1/2 — > 
+1/2 excitation and has the expectation value 

-V J dr(\u 1/2 \ 2 + \v 1/2 \ 2 ) 2 . (5.26) 

Now using that fact that u and v have spatial extent £ ~ Hvf/A and the simple BCS expres- 
sion 1 = VA^ arcsinh(c<j c /A) we arrive at the above estimate. The main conclusion is that, 
in the vortex core, it is possible for the residual interaction to give negative corrections of 
the same order of magnitude as the 'bare' excitation energy. In charged superconductors the 
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collective mode is 'plasmonized' by the coulomb repulsion. That is, there is a large positive 
correction to the energy. However this probably does not happen to the quasiparticle- 
quasihole excitation in a vortex because the corresponding translation of the vortex should 
not create a charge density fluctuation. 

In this paper we follow an alternative route and explicitly allow the order parameter to 
be dynamic, self-consistent and to affect the time evolution of the quasiparticles. It is, in 
a sense, a conserving approximation because it respects the translational invariance of the 
problem. 

To begin the calculation of conductivity we enumerate the sources of electric field and 
current density and consider only the uniform q = component. The electric field due to 
the time derivative of the background supercurrent is £ = {iujm/e)\s- There is also the 
electric field due to transverse vortex motion, N v (h/2e)zxv L , where N v is the areal vortex 
density. The trivial piece of the current density is simply the supercurrent n s e\s- Then 
there is the q = component of the current density due to vortex motion, from Eq. ( |5.20| ), 

eh 2 k 2 F± (v L - v s ) 



J g =0 — N v - 



(5.27) 



2m(A 2 /E F ) 

We have added the term — vs so that the total current is n s evs when \ F = V5. The next 
step is to write both the total electric field and current density in terms of V5 using the 
vortex equation of motion Eq. (|5.21|) . Then the conductivity tensor, a a b, may be found by 
comparing those two expressions and using the definition J a = <J a b£b- The final result for 
the conductivity tensor is rather complicated. To simplify it we again take n s = k F± /27r, 
Ep = h 2 k F± /2m and define $ = N v h/2mQo ($ is roughly H/H C 2). The longitudinal and 
Hall conductivities are 



and 
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where D = (1 + itor) 2 + V&t 2 and 



ne 



o 



^ 2 + $Ml 



1 + iur\- 
D J 



-i 



(5.30) 



We are interested in the term in the conductivity proportional to $ or the number of 
vortices when $ is small. However some care has to be taken with that limit because in 
the low frequency limit the resistivity is proportional to <3>. Now, if we throw away terms of 
order <3> 2 we have 



C.x.x ~ 



ne 2 l-$±±f^ 
mQ % + ^$ 



(5.31) 



Note that in the u — > limit this has the correct value $ 1 ne 2 r/m. In order to bring out 
the structure, we take the r — > oo limit and obtain 

-i 



ne 



imuo 



1 - 



itu 2 $ 



1 - 



ifi 2 $ 



a;T(f2o — cu 2 ) 



(5.32) 



At 'resonance', to ~ f^o, the longitudinal conductivity becomes, in this limit, completely 
imaginary. This agrees with the results of the previous subsection. 

Similarly we may look analytically at \m.a xy in the small $ limit. Imcr^y gives the 
polarization dependent dissipation. In this limit it has the value 

2ne 2 $f2 r 



Imcr^ - 



muo 



1-uj 2 t 2 + Q 2 r 2 ) +(4 



LOT) 



(5.33) 



This expression has a pole-like feature at u w Qq but it is weak because its value at the 
maximum goes as 1/r and therefore doesn't sharpen as r increases. 

In Fig. (0) we plot the real and imaginary parts of the longitudinal and Hall conductivi- 
ties as a function of frequency without the above approximations. They have been multiplied 
by uj or u 2 in order to bring out their behaviour near uj ~ Q . In general the conductivity 
at the characteristic frequency does not show dramatic behaviour. What is important is 
that Kea xx does not show the expected peak from a naive dipole resonance and lma xy does 
not show a strong peak as expected from a polarization sensitive dipole resonance. There 
is structure at low frequencies which depends on the value of $ but is not important to 
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the present discussion. It is related to the problem with taking the low frequency limit. At 
some point the resistivity becomes proportional to <3> instead of the conductivity. It is also 
the point at which the electric field due to transverse vortex motion is comparable to that 
of the time derivative of the gauge field. 



VI. PINNING 



As mentioned in the introduction, in most practical cases pinning comes into play. In this 
section we introduce the effects of pinning phenomenologically and show that a dipole-like 
resonance is recovered. First let us return to the homogeneous solution of the equation of 



motion Eq. (|5.17| ). It should be possible to cause these modes to absorb energy, however 
the equation of motion we derived apparently does not allow it to happen resonantly. It 
would seem that one requires a perturbation to 'bump' the system into one of these excited 
states. That will be the purpose of the pinning center. 

The main feature of pinning sites or bending of vortices is that there is a non-translation- 
invariant restoring force present. Thus we propose to study the effect of modifying, in a 
very general way, the equation of motion, Eq. (|5.16| ) into 



vl = v 5 + fi (v L -v s )xz v L - a 2 r . 

T 



(6.1) 



Here r is the position of the vortex and a is the characteristic frequency of a harmonic well 
we have introduced at the origin. 

Solving this equation (and keeping only terms to order 1/r) we find the steady state 
solution to be 



V L x ^ 


= iuo 
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(6.2) 



(6.3) 
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Let us now also take a « VLq. In that case the quantity in square brackets gives a resonance 
at uj w Q + a 2 /Qq. The term most important at resonance comes from the imaginary part 
of the square bracket. Let us evaluate Eq. ( |6.3| ) close to resonance by setting uj = fl Q in 
other parts of the expression to obtain 

1 



VLx 



a 2 r/(2n ) 
l + (u- n ) 2 r 2 



V 



vs. 

\ v Sy J 



(6.4) 



\ V Ly J 

The complete expression for the conductivity is messy and perhaps not worth writing 
down for this simple model. However in the spirit of subsection |V B| we may look at the dis- 
sipation in a single vortex. Setting vs y = vs x exp iO = v s j \f2 we may repeat the calculations 
leading to Eq. ( |5.22|) and Eq. fl5.23|) and obtain 



Re{v* Lx vs y ~ vl y vs x ) = f(l + sm6)vl 



(6.5) 



and 



where 



Re(v* Ly iujvsy + v* Lx iuv Sx ) = + sm9)v 2 c 



f 



a 2 r/(2n c 



(6.6) 



(6.7) 



i + (uj - n ) 2 r 2 

gives the resonance structure which corresponds to the usual dipole resonance. Note that f is 
proportional to a 2 . The dipole absorption is proportional to the strength of the translation- 
invariance-breaking perturbation. This perturbation increases somewhat the resonant fre- 
quency but need not dominate. The resonance is also fully polarization dependent and turns 
off when sin^ = —1, or, at the correct choice of handedness of the circular polarization. 



VII. DISCUSSION 



One of the main conclusions of this work, that in a pure system the single-quasiparticle 
properties of a vortex are invisible to a long wavelenth probe, is very reminiscent of the result 
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of Mattis and Bardeen |25[ regarding the long wavelength electromagnetic response near the 
bulk gap frequency. For vortices, if r — > oo we have = and there is no dissipation. 
Mattis and Bardeen found that, in the bulk, the conductivity in the long wavelenth limit 
a(q, 10 ~ 2A) — > as q —>■ unless there are impurities present to allow the violation of 
momentum conservation. There is a basic picture which is shared by these two situations. A 
quasi-particle-quasi-hole excitation can look like a translation of the condensate as a whole 
and therefore the system responds differently from what one expects. 

In this paper we have calculated various quantities without the need to choose an ap- 
propriate 'size' for the core. In previous works HQ this was required in order to fix the 
coefficient of the 'Lorentz' or 'magnus' term. The core size was chosen so that at H = H c2 
the core volume would be the total volume of the system. Here there is no particular reason 
to extrapolate to the superconducting phase boundary. 

The low temperature limit studied here may present some simplifications. A good defi- 
nition of the 'core' of the vortex would be those states closest to zero energy which move at 
velocity v^. Some intermediate states would have to be adjusted some other way to match 
the bulk velocity vg and conserve current. They will be states further away from the chemi- 
cal potential and at low temperature perhaps not important for calculating dissipation. This 
point is discussed further in the appendix. It is not clear at present exactly how crucial it 
is. The equation of motion of Bardeen-Stephen is derived assuming that the full transport 
current v# flows through the vortex core and this core is necessarily defined by some surface 
or boundary region. From the low temperature microscopic point of view the full transport 
current does not flow through the core. That is because it is unfavourable to pay the discrete 
amount of energy required to make charge fluctuations in the core. Without these charge 
fluctuations the paramagnetic current in the core is determined by the velocity of the core 
because there aren't any other states available in the vortex core. At higher temperatures, 
it is favourable to allow charge fluctuations in the core and thus it is possible to increase the 
current flowing there. 

Future work should concentrate on detailing the microscopic theory of pinning and vortex 
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bending since that would be most important for the experiments. It would be very nice if the 
de-pinning frequency correlation presented in table | could be explained. It would also be 
useful to work out the quasiparticle structure of a vortex in a d-wave superconductor That 
involves a non-local BdG equation. One could also check how the discreteness of the lattice 
affects the microscopic structure. In general, calculations of the microscopic strcuture will 
be sensitive to the nature of the pairing. Thus knowledge of the microscopic strcuture of 
vortices may give us some insight into the mechanism of high temperature superconductivty. 

VIII. CONCLUSION 

In this paper we have studied the microscopic dynamics of vortices as a function of 
frequency. We have derived an equation of motion and calculated the electrodynamical re- 
sponse at q — 0. For free vortices we find that the single-particle character is essentially 
invisible unless translation invariance is broken. Thus the clear, albeit broad, response in 
optical transmission of Karrai et al. is perhaps due to pinned vortices. Another interesting 
possibility is that interaction with the discrete lattice is responsible. The possibly d-wave 
nature of the paired state in high temperature superconductors should not affect this con- 
clusion. We believe that the large width observed experimentally may come from intrinsic 
scattering from whatever breaks translation invariance. 
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MATRIX ELEMENTS 

Throughout this paper we use matrix elements of the quasiparticle wavefunctions. They 
are difficult to calculate because the detailed behaviour of the wavefunctions depends non- 
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trivially on the radial dependence of |A(r)| which in turn is determined self-consistently. In 
this appendix we derive approximate matrix elements of r and V between the low energy 
quasiparticle states. The approximation we make is that the main contribution to these 
matrix elements comes from close to the center of the vortex where the functions u(r) 
and v(r) can be approximated by Bessel functions and the order parameter A(r) is small. 
It is assumed that matrix elements of the low energy states are not very sensitive to the 
exact behaviour of the wavefunctions near the core boundary because the functions u and 
v are suppressed exponentially when the magnitude |A(r)| becomes substantial. Boundary 
conditions have entered implicitly through the energies of the quasiparticle levels. 
We shall make use of the Bessel function identities 

J v -i(z) + J v+1 (z) = 2 fJ„(z), 

J v - 1 (z)-J t + 1 (z) = 2±J v (z). (1) 

and suppose that u(r) and v(r) satisfy similar identities at least where they contribute most 
to the matrix element. To be specific, in the calculation of matrix elements of V we use 

^u, = ^ (e*V-i " e-% + i) 
~« M = - Ht ( e *Vi + e"*Vn) (2) 

and exactly the same equations with u replaced by v. 6 is the polar angle. The matrix 
elements of 

n ■ V = cos(fl - 6 ) — - sin(# - 6 )— (3) 

where n is a unit vector pointing at angle 9 to the x-axis, are thus 

d 3 r^ ±1 n • VVv = ^e* + W 2 » (4) 

and zero for states differing by other angular momenta. The lack of dependence on /i is a 
result of the approximations we have made. 
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Another approximation we have made is to only consider matrix elements between low 
energy states. We consider in this paper mostly operators that contain angular momentum 
I = ±1 and therefore amongst the low energy states it is clear which matrix elements are 
non-zero. The approximation made in Eq. (0) does not affect that. There are, however, 
high energy scattering states (of energy A or higher) which have non-zero matrix elements. 
These have been ignored. It is not clear whether this is a safe approximation or not. 

Now let us consider the matrix elements of a z hr. To do this we consider the commutator 
[H, a z h ■ r] where H is the BdG Hamiltonian with A = 0. We have 



J V^'n ■ xi\) v = J ipl[H, a z h ■ r]^ 



711 



-n • p + 2n 



-A 
A* 



(5) 



The second term in the square brackets is a smaller contribution because whenever u and 
v are substantial, |A| is small so in our approximation we shall ignore it. We are then left 
with the matrix element of the momentum operator which we have just evaluated above. 
The matrix element is 



-hk 



2m (ep - ( 



fJ,,U—l c VfI,,U+l 



(6) 



In section |V] we encounter the matrix element 



1 



5A* 



(7) 



/ 



where 5A = —Sr' ■ VA. It may be found by utilizing translation invariance if fi ^ v. 
Let Sip u (r) = —Sr' ■ Vip u (r) be the change in the eigenfunctions upon displacing them by 
an amount 5r' Q . Then substituting 5ip v and SA into the Schrodinger equation, Eq. (|2 . 2|) , 
subtracting the same equation with 5r' = 0, keeping terms to first order in 8r' , multiplying 
on the left by Vv( r ) an< ^ integrating by parts results in 



W, 



(8) 
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This can be further simplified making the assumption again that u(r) and v(r) resemble 
Bessel functions in the relevant spatial region yielding 

5rk F ± 



5v' ■ VMr'o) 



2i 



e^o+f )^_ l(r ) + e -Wo+5 )^ +1 ( r ) 



(9) 



where 5r makes an angle with the x axis. In terms of the displacement 5r' Q the matrix 
element is, 



W vv = 



Sr'k 



0^F± 



2i 



(10) 
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FIGURES 

FIG. 1. (a) Dissipation (normalized to the high frequency limit) as a function of frequency 
(normalized to Oo) when £Iqt = 1. Solid line: circular polarization of vs; dashed line: opposite 
sense of circular polarization, (b) The same quantities but for Qot = 3. 

FIG. 2. (a) Real part of the longitudinal conductivity multiplied by frequency oj (for clarity), 
Reu)a xx , as a function of frequency (normalized to Qq). Solid line: <E> = 0.1; dashed line: $ = 0.2 . 
(b) Same plots for the imaginary part multiplied by lu 2 for clarity, Imio 2 a xx . (c) Real part of the 
Hall conductivity, Reoj 2 a xy . (d) Imaginary part, lmuj 2 a xy . 
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TABLES 

TABLE I. Observed crossover frequencies for vortex de-pinning and quasiparticle energy level 
separations in vortex cores (A 2 /Ep, where A is the bulk gap and Ep is the Fermi energy) for 
various samples of low temperature superconductors. 

T C {K) Observed frequency (MHz) Energy level separation (MHz) 

Pbln g 1.7 4 2 

Pbln g 1.7 5 2 

NbTa g 4.2 15 20 

Pbln @ 1.7 7 2 

Pbln @ 1.7 8 2 

NbTa Q 4.2 26 20 

NbSe 2 pi] 7.2 100 1000 



33 



